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AN ANALYSIS OF THE DAHL FRICTION MODEL 
AND ITS EFFECT ON A CMG GIMBAL 
RATE CONTROLLER 

I. INTRODUCTION 


Rocent design work on the fine pointing control system for the Large 
Space Telescope (LST) has generated a keen interest in the low level torquing 
characteristics of momentum exchange devices, particularly medium si/.e con- 
trol moment gyros (CMGs). This attention to CMG actuators stems from the 
desire to find a single actuator that is adequate for all pointing and maneuver 
modes of the LST. The definition of low level torquing for the purposes herein 
is that level of torque resolution required to point the LST to ~ 0.001 arc sec, 
or approximately 0.01 N-m for a control system natural frequency of 1 11 k. The 
difficulty in achieving tills torque resolution becomes apparent when one calcu- 
lates the associated gtmbal rates. For a 200 N-m-s CMG, as an example, the 
gtmbal must be accurately processed at a rate ol 5 x 10 ■' rad/s, an angular rate 
nearly an order of magnitude lower than that of the hour hand ol a clock. Hie 
problems in building a CMG and its gimbal control system for good low signal 
performance are primai’ily caused by frictional torques in the gimbal pivots. 
Frictional torques in bearings and brushes generally exhibit a nonlinear charac- 
teristic at small angular velocities and, consequently, complicate gimbal rate 
control for small rates about zero. 


lire precise nature of the frictional mechanism is the subject of some 
controversy. Since applications that required such extremely low torque levels 
were nonexistent until recently, only a small amount of Information on the 


subject, either theoretical or experimental, has been published. Consequently, 
the classical stietion characteristic was used in the initial modeling work. ‘Hits 
characteristic, of course, made the medium size CMG undesirable as an 
actuator for a fine pointing system, and it was not long before the validity of the 
stietion model began lobe questioned. Based on Investigations bv Bald 4 and 
some recent experimental work at Sperry Flight Systems in Phoenix, Arizona, 


♦Bald, P. R.J A Solid Friction Model. Aerospace Report No. TOR-birrs 
(:tl07-lH)-l, May IflKH. 


tt Ih now generally felt that the frictional torque In a glmbnl pivot Is not at all 
characterized by the classical sttctlon model but Is more closely represented 
by a nonlinear spring, with torque botng a function of the glmbnl displacement, 
In his roport Dahl discusses the various mechanisms for frictional forces 
arising from various surface conditions In both rolling and rubbing eontnet. 
Special emphasis Is glvon to the modol for friction In ball bearings, Including 
experimental data which show close agreement with the analytical model, 

Whilo this agree mont In ttsolf lends strong support to the validity of the model, 
Dahl’ s representation also soems to Incorporate nil the phenomena that one 
might intuitively consider essential to the friction mechanism, 

The objective of this document is to extend the interpretation of what the 
Dahl model does; nothing new is added to the basic theory. There is n section 
that describes the model and another one, entitled Analytical Results, wherein 
some interesting dynamical characteristics become apparent as the result of 
equation manipulations. Another section gives the results of a digital simula- 
tion of the equations of motion for rotating mass supported by n roller bearing 
pivot and influenced only by frictional torque. The Inst section indicates how 
frictional torques affect the design of a glmbnl rate control system. 


II. DESCRIPTION OF THE SQUARE LAW FRICTION MODEL 

If Tj(0) is defined as the frictional torque, its time derivative can be 
written as 

dt (10 dt * ' ^ 


where 0 is the shaft angular displacement, IT, then, one can define the function 
dTj/dfl for a given bearing pivot, the frictional torque can be computed. It is 

apparent, then, that the choice of dT^/dfl determines the model, Dahl has 

postulated the following form for this function: 



*> 


f 


»■* 


(ho ronfumnblonoflfl of which In supported by experimental data, A block diagram 
representing equations (l) and (l i) Is shown In Figure l « Is Uie maximum 

frlctlonnl torque that enn occur nnd In equivalent to what In commonly termed 
running friction, 



To show the operation of the model, one shoi.d determine what happens 
within the block diagram as the input, 0, is varied. It is important to note, 
first of all, that dT^/dfl , the output of the square law block, is always positive. 

Assume that initially all variables are zero, then 0 goes positive. Since, for 
T f = 0, dTj/dO = yT* 0 , T f will have a finite positive value nnd T ( . will begin 

to increase. Because ‘ is positive the switch is in the positive position; and 
as T increases, cITj/dO becomes smaller. This process continues until 

T = T f() when dT^dfl = 0 and T - 0 so that T f can no longer increase. If 

0 becomes negative, the input to the square law block changes sign and, at the 
instant 0 goes negative, dT^/d 0 jumps from zero to tr T j* n • Huis, T f Is 

negative and the integrator begins integrating down from . It is inU*restlng 


to note at this point that the frictional torque is in a direction to increase the 
angular velocity in the negative direction; or colloquially, the friction torque 
is aiding the motion. This kind of phenomenom is generally not associated with 
frictional torques; and, if one defines friction as an energy loss mechanism, 
then it is obvious that the Dahl model incorporates more than energy loss in the 
pivot. The integrator then continues to integrate from T to -T^ . 





In his report Dahl Indicates that Tj, has Home randomness nHHoelnted 

with It due to the nonuniformity of the two contacting HurfneeH and he has Hhovvn 
how thtH might be incorporated Into the model, This port of the model, how- 
ever, Is not considered In this report. 


III. ANALYTICAL CONSIDERATIONS 


In the case of a mass free to turn In a pivot and Influenced by frictional 
torques, the equations defining its motion ares 



and 



(3) 


(•0 


where I Is the moment of inertia of the rotating mass about the pivot axis,' 

* 

Now, consider small motions about the equilibrium point, T = 0 ~ 0, so 
that the equations can be linearized about this point. If the square function on 
the right-hand side of equation (4) is expanded in a Taylor series, all terms 
except the first are cither zero or second order, and equation (l) reduces to 


V 


yT io 0 
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which is directly tntograble. hints, equation (a) becomes 


\0 + yT‘ () 0 


(•'>) 


the oquntlon for n linear oscillator with y'P* being the spring rate and natural 
frequency to given an 


to 


s 



1 


00 


Hence, to first order, the friction model behaves ns an oscillatory second order 
system with no damping. The absence of a damping term In the linearized 
equations doos not, howover, preclude the possibility of asymptotic stability. 

It Is also of Interest to compute T f ns a function of 0 and n where 
n a A sin to 0 t. Kq nation (4) can then be Integrated In the form 


T f dT f 
! VtT f -T f0 C 


t 

f A sin to t dt 
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(?) 


for 0 ^ to Q t =s 7T. Integration of equation (7) gives 
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Hie proper value of T. for a given A and to n can bo found from equation (s) 

l . 

from tho condition that T = T„ Tor to 0 t = ff , and the ratio T /T has the 

1 l 'j 


f 


form 







( 0 ) 


Kquntton (8) 1ms boon plotted as n function of 0 and 0 In Figures 2 and U. 
The angular velocity, o*, Is a sinusoidal function with a constant amplitude 
A = 0,002 rad/s and a> 0 Is varied as a parameter. 


These curves, plotted from simulation results, appear in Dahl* s report, 
and he makes the observation that ns co 0 decreases, the shapes of the Tj. 

versus 0 curves approach the form of the classical stlctton characteristic. 
Tills can be seen from equation (8), since ns co 0 -* 0, T«/T« 0 -* 1 for an 

arbitrary value of t. On the other hand, as to 0 Increases, the curves become 
more elliptical in shape. Tills Is Indicative of die fact, for T^« 1, T^ Is a 

linear function of 0 as shown In equations (4) and (5) which will produce nn 
ellipse. Likewise die curves of T^ versus 0 approach a straight line through 
the origin. 


The quantify, 1/yT^ , has the unit of rad and can be Interpreted as a 

characteristic angle for the friction model. The relative size of A/x' 0 , the 
amplitude of tile angular motion, widi respect to the characteristic angle deter- 
mines die shape of the curves In Figures 2 and 3, It can be shown, for Instance, 
that the two Inflection points In die graphs for small a> 0 ' s disappear for 
T y (A/'jo 0 ) e 2. Hence, die shapes of the curves change most significantly 

for amplitudes that are approximately twice the characteristic angle. 


IV. SIMULATION RESULTS 


In order to verify the analytical results of die previous section, a digital 
simulation was developed for the gimbnl control system depleted In Figure I. 
'Hie equations of motion for this system arc the followings 


(5 






Figure 4. Gtmbnl rate controller. 


10 = K (o\ - 0 ) +. K (0 -0) - T 
p i II t 


( 10 ) 


Vi - WV 


and 



00 


( I ~ J 


wlicrc l is the effective Inertia of the gtmbnl, K and are control system 
gains, and o ^ Is the commanded gtmbnl rate, 1 

Since Initially Interest lies In the dynamics due to frictional torques 

alone, K anil K. are set equal to zero. If, then, the rotor Is given an Initial 

pi 

velocity, It should come to rest under the Influence of the frictional torque. 








rho motion during thlfl time will Indicate the general characteristics of tlu> 
friction model, Figure 5 shows results for nn Initial n ~ n.nOM rnd/s, and 
V N~m, n,1f l T r ’ d, Hfi v | o'* (N- in rnd)"*j the mo values of T^ nnd 

T were men sured experl mentnlly by Dahl. The Inertia, I, was chosen ns r, 
.slug ft’ which Ih typical of the effective glmbnl Inertia of n medium .size CMC), 
Figure 5 shows Hint for the first quarter of a second, until 0 ehnnges sign, 

'r f T f() . The motion thereafter Is similar to that of an underdamped second 

order system, 'Hie period of the oscillatory motion taken from Figure 5 Is 
13. « seconds while that computed from equation ((’.) Is 13.05 seconds. Thus, the 
spring mechanism Is evident. It Is also noteworthy that as the amplitudes 
become smaller, the damping decreases, as was found above. 


Figure (1 shows similar results for a ease with T = 0. 27 N-m and 

y~ 1.48 x io r » (N-m rad) -1 . These values are the result of fitting some 
experimental data taken on a CMG pivot. Those numbers have Increased the 
natural frequency considerably with the measured period equal to (). Ills second 
and that calculated from oquution (4) 0. 135 second. If one compares Hie plots 
of 0 and I in 4 tgures 5 and 0, it appears that the two variables are in phase 


and that T f is sinusoidal in nature. However, a closer look at T subsequent 

to the times when 0 is changing sign Indicates that T decreases more rapidly 

than it would If it wore a true sine wave. This Is duo to the switching from a 
low to a higher gain in the friction model and has the effect of producing damp- 
ing. As the amplitude decreases, this effect is seen to diminish, since the 
change in gain at the switching instant becomes loss and loss. 


The characteristics shown in Figures 2 and 3 wore verified by command- 
ing a sinusoidal input for n ( through the control system with K - kv 1 nnd 

0 I 

Kp = -M. The results compare well with those calculated in the previous 
section. 


V. THE EFFECT OF THE FRICTION MODEL ON THE 
GIMBAL RATE CONTROL SYSTEM 

If glmbnl pivot friction is neglected, the glmbnl rale controller Is given 
by equations ( 10) and (It) with T ( . - (). Hie natural frequency, :o _ , and damp- 
ing laetor, f, , for (lie controller are give’ll by 
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Motion caused by frictional torques 






Motion caused by frictional torques. 
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While the design natural frequency of the control loop depends on the particular 
application, generally it is above 5 Ilz. In Section HI it was shown that there 
was a natural frequency associated with the pivot friction loop given by the 
equation 




I 


For the values of y and T fQ in Dahl’ s report, co f = 0. 28 Hz as shown in 

Figure 5. Figure 6, on the other hand, shows that, for the values obtained 
from a curve fit to CMG data, 7.4 Hz. Hence, it is conceivable that 

could be similar in magnitude to u >^ . To what extent, then, is such a situation 

detrimental to the system 1 s performance? 

Figures 7 through 10 show the step response for the system defined by 
equations (lo), (ll), and (12) for various parameter values. The step com- 
manded gtmbal rates, 0 ^ of r,x 10"' 1 rad/s or 5* lo- r > rad/s are used to show 

the effect of commanded magnitude on the system response. Figure 11 shows 
that, for the values of y and T^ in Dahl* s report, the natural frequency of 

the friction loop, , and the frictional torque are so small that even for the 

low level command the friction has very little effect} the response with and 
without friction are very nearly identical. Figures 8 and 9 show that for 
T f0~ °*“ 7 N “ m y ~ t,,l8x 1° B f N-m rad | the response is affected significantly 

by the magnitude of 0^ in the range of interest. In the responses of Figures Hi 

and 1 1 the eontiol system natural frequency has been increased by a factor of 
2, and still there is serious deviation front the desired response. 
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Figure 7. Gimbal rate controller step response. 
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Figure 9. Gimbai rate control 




Figure 10. Gimbal rate controller step response. 
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Figure 11. Gimbul rate controller stop response. 






In Figures 12 through \<\ 0^ in n sinusoidal function with n frequency of 
2n rnd/s , The loss In gnln In apparent ns the amplitude In decreased. 


VI. SUMMARY AND CONCLUSIONS 


The basic characteristics of the Dahl friction model have been shown by 
means of analytical and simulation methods. The effects of friction, repre- 
sented by the Dahl model, on a CMC! gimbal rate control system was investigated 
by digital simulation. The conclusion from these simulation results is that 
gimbal pivot friction can have a significant effect on the gimbal rate control 
system, The magnitude of the problem this presents depends on the character- 
istics of the actual pivot. It would appear from this preliminary look that one 
solution is to insure that the control system natural frequency is higher by some 
prescribed amount than the natural frequency of the friction loop. 
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